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Through the measurements of magnetic field dependence of specific heat in La2- x Sr x CuOi in 
zero temperature limit, we determined the nodal slope % of the quasiparticle gap. It is found that 
wa has a very similar doping dependence of the pseudogap temperature T* or value A p . Mean- 
while the virtual maximum gap at (ir, 0) derived from i>a is found to follow the simple relation 
A q = 0.46/cbT* upon changing the doping concentration. This strongly suggests a close relation- 
ship between the pseudogap and superconductivity. It is further found that the superconducting 
transition temperature is determined by both the residual density of states of the pseudogap phase 
and the nodal gap slope in the zero temperature limit, namely, T c w /?VA7n(0), where 7n(0) is the 
extracted zero temperature value of the normal state specific heat coefficient which is proportional 
to the size of the residual Fermi arc k arc . This manifests that the superconductivity may be formed 
by forming a new gap on the Fermi arcs near nodes below T c . These observations mimic the key 
predictions of the SU(2) slave boson theory based on the general resonating-valence-bond (RVB) 
picture. 

PACS numbers: 74.20.Rp, 74.25. Dw, 74.25. Fy, 74.72.Dn 



INTRODUCTION 

Since the discovery of the cuprate superconductors, 19 
years have elapsed without a consensus about its mecha- 
nism. Many exotic features beyond the Bardeen-Cooper- 
Schrieffer theory have been observed. One of them is the 
observation of a pseudogap in the electron spectral func- 
tion near the antinodal points (tt,0) and (0,7r) at a tem- 
perature T* >> T c 1]. In a conventional BCS supercon- 
ductor, this gapping process occurs simultaneously with 
the superconductivity at T c . It has been heavily debated 
about the relationship between the pseudogap and the su- 
perconductivity in cuprates. One scenario assumes that 
the pseudogap A p marks only a competing or coexisting 
order with the superconductivity and it has nothing to do 
with the pairing origin. However another picture, namely 
the Anderson's resonating-valence-bond (RVB)[2( model 
(and its offspring) |3( predict that the spin-singlet pair- 
ing in the RVB state (which causes the formation of the 
pseudogap) may lend its pairing strength to the mobile 
electrons and make them to naturally pair and then to 
condense at T c . According to this picture there should 
be a close relationship between the pseudogap and the 
superconductivity. 

In order to check whether this basic idea is correct, we 
need to collect the information for the pseudogap and the 
superconducting energy scale, especially their doping de- 
pendence. The pseudogap values A p (or its correspond- 
ing temperature fc^T* ~ A p ) and its doping dependence 
have been measured through experiments 1]. To deter- 
mine the superconducting energy scale, we note that the 
normal state Fermi surface is formed by four small arcs 



near the nodal points |4|. As temperature is lowered be- 
low T c , a new gap opens on these arcs. To illustrate this 
point more clearly, in Fig.l we present a schematic plot 
for different gaps or energy scales. The dotted line repre- 
sents the gap structure of the normal state, assuming the 
presence of Fermi arcs near the nodal points. The region 
of zero gap corresponds to the Fermi arc. The dash line 
and the solid line represent two possible gap structures 
for superconducting state at T = 0. The solid line is 
the standard d-wave gap with maximum gap value A p at 
("7T, 0) and (0, 7f). From this picture, we see that the nodal 
gap slope, which is defined as «a = [dA a /d6] no( i e /hkF, 
can be used to determine the superconducting energy 
scale. The relationship between the nodal gap slope va 
and the maximum pseudogap A p remains to be a big puz- 
zle. In particular, the two quantities may be independent 
of each other if the superconductivity is not induced by 
the formation of the pseudogap. Therefore to measure 
the nodal gap slope near nodal point in the zero temper- 
ature limit becomes highly desired. When combined with 
the known results on the pseudogap A p , this will allow 
us to detect the relation between the pseudogap and the 
superconductivity. 

Some previous results using, for example, angle- 
resolved photo-emission (ARPES) 5] or superfluid den- 
sity seem to be inconclusive due to either energy resolu- 
tion (ARPES above lOmeV) or unexpected difficulty in 
analyzing the data (e.g., a so-called Fermi liquid correc- 
tion factor (Xfl is inevitably involved in analyzing the 
low temperature data of superfluid density). In this pa- 
per, we report the evidence of a proportionality between 
the nodal gap slope «a and the pseudogap temperature 
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FIG. 1: Schematic plot for the pseudogap energy, supercon- 
ducting energy scale and nodal gap slope. The solid line rep- 
resents a standard d-wave gap A = A p cos26 with maximum 
gap value A p near (ir,0) and 6 the angle starting from k x . 
The dotted line shows the pseudogap near nodes if the su- 
perconductivity would be suppressed completely ( based on 
the Fermi arc picture ). The dashed line shows a possible 
quasiparticle gap near nodes. The nodal gap slope is defined 
as ha = [dA 3 /d8] no d e /hkF. The nodal gap slope va and the 
maximum gap A p near (n, 0) may not be related if the super- 
conductivity (which controls the gap structure near nodes) 
has nothing to do with the pseudogap. 

T*. Remarkably a simple relation, namely A q — 0.46T*, 
between the virtual maximum quasiparticle gap (A q ) 
derived from va and the pseudogap temperature T* is 
found. We also find that T c is determined by both the 
nodal gap slope va and the size of the Fermi arcs (k arc ) 
in the underdoped normal state. Both observations are 
anticipated by the SU(2) slave boson theory based on 
the general RVB picture. 



EXPERIMENT 

We determine the properties of the nodal quasiparti- 
cles by measuring the low temperature electronic specific 
heat. The La2-xSr x CuO± single crystals measured in 
this work were prepared by travelling solvent floating- 
zone technique. Samples with seven different doping 
concentrations p=0.063(T c =9K, nominal x=0.063, post- 
annealed in Ar gas at 800°C for 48 hrs ), 0.069(T C =12K, 
as-grown sample with x=0.063), 0.075 (T C =15.6K, nom- 
inal x=0.07 and post-annealed in O2 gas at 750°C 
for 12 hrs), 0.09 (T C =24.4K, as grown, x=0.09), 0.11 
(T C =29.3K, as grown, x=0.11), 0.15 (T C =36.1K, nomi- 
nal x=0.15), 0.22 (T C =27.4K, nominal x=0.22) have been 
investigated. The quality of our samples has been char- 
acterized by x-ray diffraction, and R(T) data showing 
a narrow transition AT C < 2 K. The samples have also 
been checked by AC and DC magnetization showing also 
quite narrow transitions. The full squares in Fig. 6 rep- 



FIG. 2: Field dependence of A7 = \C{H) - C*(0)]/T normal- 
ized by the data at about 12 T in zero temperature limit. It 
is clear that Volovik's y/H relation describes the data rather 
well for all samples. This indicates a robust d — wave super- 
conductivity in all doping regimes. 



resent the transition temperatures of our samples. The 
heat capacity presented here was measured with the re- 
laxation method based on an Oxford cryogenic system 
Maglab-EXA-12. In all measurements the magnetic field 
was applied parallel to c-axis. As also observed by other 
groups for La— 214 system, the anomalous upturn of C JT 
due to the Schottky anomaly of free spins is very weak. 
This avoids the complexity in the data analysis. Details 
about the sample characterization, the specific heat mea- 
surement, the residual linear term and extensive analysis 
are reported in a recent paper 0. 

It has been widely perceived that the pairing symme- 
try in the hole doped cuprate superconductors is of d- 
wave with line nodes in the gap function. In the mixed 
state, due to the presence of vortices, Volovik || pointed 
out that supercurrents around a vortex core lead to a 
Doppler shift to the quasi-particle excitation spectrum. 
This will dominate the low energy excitation and the 
specific heat (per mol) behaves as0, C vo i = AVH 
with A cx 1/va- This square-root relation has been 
verified by many measurements which were taken as 
evidence for d-wav e sy mmetry, for example by specific 
heat [lol ITll fl^. Il4l IT5L IT ^I . thermal conductivity [l7| . 
tunnelling (to measure the Doppler shift of the Andreev 
bound states) 0] , etc. In this way one can determine 
the nodal gap slope (ua)- Since the phonon part of the 
specific heat is independent on the magnetic field, this 
allows to remove the phonon contribution by subtracting 
the C/T at a certain field with that at zero field, one 
has A 7 = AC IT = \C{H) - C(0)]/T = C vo i/T - aT 
with a the coefficient for the quasiparticle excitations of 
a d-wave superconductor at zero field (C e = aT 2 ). In 
the zero temperature limit A7 = C vo i/T = A\/~H is an- 
ticipated. 
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FIG. 3: (a) The typical original data of A7 vs. T for the 
underdoped sample p = 0.069 at different magnetic fields, 
(b) The same set of data plotted as Aj/yH vs. T. One 
can clearly see that in zero temperature limit A^/i/H is a 
constant for all fields implying the validity of the Volovik's 
relation A7 = AyH. From here one can also determine the 
value A which is about 0.28rnJ/molK 2 T 5 as marked by the 
thick bar. 



RESULTS AND DISCUSSION 

In order to get A7 in the zero temperature limit, we 
extrapolate the low temperature data of C/T vs. T 2 (be- 
tween 2K to 4K) to zero K. The data taken in this way 
and normalized at 12 T are presented in Fig. 2. It should 
be mentioned that the similar data have been published 
in our previous paper 0- Here for clarity we present 
the data again with more detailed analysis. It is clear 
that the Volovik's y/~H relation describes the data rather 
well for all doping concentrations. This is to our surprise 
since it has been questioned whether the Volovik relation 
is still obeyed in the underdoped regime especially 
when competing orders are expected to anpear |l9ll2fll2 1| 
and impurity scattering is present. We attribute the 
success of using the Volovik relation here to three rea- 
sons: (1) We use A7 = [Cfl-iic — Ch=o\/T instead of 
using A7 = [Ctji| c — Ch± c ]/T. The latter may in- 
evitably involve the unknown DOS contributions from 
other kinds of vortices (for example, Josephson vortices) 
when H _L C, (2) The contribution from a second 
competing order to A7 may be small compared to the 



Volovik's term in the zero temperature limit. This is 
reasonable when considering a contribution to the heat 
capacity by the competing order as ~ T v with v > 1. 
For example, the specific heat due to the spin correlation 
in 2D anti- ferromagnetic phase is ~ T 2 . In zero temper- 
ature limit this term goes away. (3) The DOS induced 
by the Doppler shift effect in our experiment is much 
stronger than that induced by the impurity scattering. 
We will further address this point in the forthcoming 
discussion. To have a self-consistent check of the y/~H 
relation found in the zero temperature limit, we plot the 
raw data of Aj/y/H vs. T at finite temperatures. A 
typical example for the very underdoped one (p = 0.069) 
is shown in Fig. 3 (a) and (b). One can see that in the low 
temperature region the data Aj/y/H scale for all fields 
ranging from 1 T to 12 T, showing the nice consistency 
with the relationship A7 cx \/~H for this sample in the 
zero temperature limit. From here one can also deter- 
mine the prefactor A in A7 = A-\/H (here for example, 
A = Q.28mJ/molK 2 T - 5 for p=0.069) and then compare 
backwards to the value determined from the data shown 
in the main panel leading to of course the same value. 
The same feature appears for all other doping concentra- 
tions. For clarity they will not be shown here. 

It is clear that the Volovik's \f~H relation describes the 
data rather well for all doping concentrations. This suc- 
cessful scaling of A7 vs. yH makes it possible to derive 
the pre-factor A, and one can further determine the nodal 
gap slope va- Fig. 4(a) shows the doping dependence of 
the pre-factor A. The error bar is obtained by fitting the 
extracted zero temperature data to A7 = A\/~H. For a 
typical d-wave superconductor, by calculating the exci- 
tation spectrum near the nodes, it was shown that [lf3 | 
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here l c = 13.28 A is the c-axis lattice constant, V mo i 
=58 cm 3 (the volume per mol), a p a dimensionless con- 
stant taking 0.5 (0.465) for a square (triangle) vortex 
lattice, n=2 (the number of Cu-0 plane in one unit cell), 
$0 the flux quanta. The va has then been calculated 
without any adjusting parameter (taking cv p =0.465) and 
shown in Fig. 4(b). It is remarkable that v& has a very 
similar doping dependence as the pseudogap tempera- 
ture T* , indicating that «a « T* a A p . If convert- 
ing the data t>A into the virtual maximum quasiparticle 
gap (A 7 )|iij via «a = 2A q /hkp, here kp ~ Tr/y/2a is 
the Fermi vector of the nodal point with a — 3.8A (the 
in-plane lattice constant), surprisingly the resultant A q 
value [shown by the filled squares in Fig. 4(b)] is related 
to T* in a simple way (A g ~ 0.46fc_eT*). It is impor- 
tant to emphasize that this result is obtained without 
any adjusting parameters. Counting the uncertainties 
in determining T* and the value of a p , this relation is 
remarkable since A q and T* are determined in totally 
different experiments. Because va (or A g ) reflect mainly 
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FIG. 4: (a) Doping dependence of the pre-factor A determined 
in present work (open circles). Here the point at p=0.19 was 
adopted from the work by Nohara et al. on a single crystal [l^l . 
(b) Doping dependence of the pseudogap temperature T* 
(open symbols) summarized in Ref.l (see Fig. 26 there) and 
our data da (solid line). T*-susceptibility refers to the pseu- 
dogap temperature determined from the maxima in the static 
susceptibility, and T* — p to the temperature at which there is 
a slope change in the DC resistivity (all in La2- x Sr x Cu04). 
Above T* — p the resistivity has a linear temperature de- 
pendence. The full squares represent the calculated virtual 
maximum quasi-particle gap A 9 derived from va without any 
adjusting parameters. Surprisingly both set of data are cor- 
related through a simple relation (A q « 0.46fcsT*) although 
they are determined in totally different experiments. This re- 
sult implies a close relationship between the pseudogap A p 
and the nodal gap slope va ■ 



the information near nodes which is predominantly con- 
tributed by the superconductivity, above discovery, i.e., 
va oc T* oc A p (or A q ~ 0A6k B T*) strongly suggests 
a close relationship between the superconductivity and 
the pseudogap. A similar conclusion was drawn in un- 
derdoped Y Ba^Cu^O,, by analyzing the low temperature 
thermal conductivity|22(. If the pseudogap is supposed 
to be caused by the formation of the RVB state[2j, our 
results here point to a fact that the RVB singlet pairing 
may be one of the unavoidable ingredients for supercon- 
ductivity. It remains to know whether this conclusion 
holds also for the electron doped samples since so far it 
is not clear yet whether the pseudogap exists in these 
N-type samples. 

In above discussion, we see the consistency between 
our low temperature specific heat data and the Volovik's 



FIG. 5: (a)The doping dependence of 7n(0) derived from spe- 
cific heat (open squares |24||') based on entropy conservation, 
ARPES (open circles^), and Knight shift (up-triangles^) 
in La,2-xSr x CuOi. The solid line is a fit to the data with 
7n = C(P -PcT yielding C = 182.6, p c = 0.03, r] = 1.54. (b) 
A comparison between karc calculated from 7«(0) and 2n/a. 
One can see that the k arc becomes smaller than 2n/a in the 
underdoped region showing the relevance of the Fermi arcs in 
the pseudogap phase. 



square root relation A7 = As/H. This seems surpris- 
ing since the temperature range we considered here is 
about several Kelvin. At such an energy scale, the 
impurity scattering will strongly alter the DOS in the 
low energy region by generating some new quasiparti- 
cles. However, by applying a magnetic field, the Doppler 
shift of the quasiparticle excitation spectrum will con- 
tribute a new part to DOS. This energy shift is actu- 
ally not small comparing to the temperature. We can 
give a simple estimation on the energy shift AE. It is 

known that AE = ckfl^/^t^, here Ib is the mag- 
netic length which is defined as Ib = (ch/eB) 1 / 2 , and 
0<a^L<lisa Fermi liquid correction term. Taking 
vp = 2.73x10 7 cto/s[23|. we have AE 



3.67a FLy /B/lT 
meV. For example, taking the maximum field (12 T) in 
our experiment, we get AE — Ylflap^meV which is ac- 
tually a relatively big energy scale compared to the tem- 
perature T since cvfl ~ 1< This may explain why the 
Volovik's simple square-root relation A7 = As/H can 
be easily observed in our single crystals with inevitable 
certain amount of impurities. 

In the following we will investigate what determines 
T c . Bearing the doping dependence of in mind, it is 
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easy to understand that v^TiUf should not be a good es- 
timate of the superconducting energy scale for the under- 
doped samples since the T c and va have opposite doping 
dependence. The basic reason is that the normal-state 
Fermi surfaces are small arcs of length k arc near the nodal 
points. The superconducting transition occurs by form- 
ing extra gaps on the Fermi arcs. So the effective su- 
perconducting energy scale should be estimated as E s ~ 
7}VAhk arc . From the normal state electronic specific heat 
C e ie = ~/ n T, we have 7„ = Ank%k arc V mo i/hv F l c . Assum- 
ing E s r*> ksT c we find 



&nk 3 B V mol 



(2) 



where a s is a dimensionless constant in the order of unity, 
vf is the nodal Fermi velocity normal the Fermi sur- 
face. The value of 7™(0) can be estimated from specific 
heat0 IH, or indirectly by ARPESjH or NMR^. 
Here we take the values for 7 n (0) summarized by Mat- 
suzaki et al.|24| and fit it (in unit of mJ/molK 2 ) with a 
formula j n = ({p~p c ) n yielding ( = 182.6, p c — 0.03, rj = 
1.54. In Fig. 5 we present the doping dependence of the 
zero-temperature specific heat coefficient 7n(0) and k arc . 
One can see that k arc becomes smaller than 2ir/a in 
underdoped region showing the self-consistency of the 
picture of Fermi arcs. In Fig. 6 we present the dop- 
ing dependence of the truly measured T c (filled squares) 
and the calculated value (open squares) by eq.(2) with 
(3 = 0.7445K 3 mols/ Jm. In underdoped region, the truly 
measured and calculated T c values coincide rather well 
implying the validity of eq.(2). In the overdoped region, 
7„ will gradually become doping independent, therefore 
one expects T c ex i>a- So the energy scale of the supercon- 
ductivity is not given by VA^kp ~ A p , but by \vAhk arc 
or more precisely by eq.(2) in the underdoped region. 

To have a framework about the experimental results, 
in the following, we will review one particular expla- 
nation based on the slave-boson approach. Within the 
SU(2) slave-boson theory 0, the pseudogap metallic state 
is viewed as a doped algebraic spin liquid f ASIA |2s| . 
A doped ASL is described by spinons (neutral spin- 1/2 
Dirac fermions) and holons (spinless charge-e boson) cou- 
pled to a U(l) gauge field. Due to the attraction between 
the spinons and the holons caused by the U(l) gauge 
field, a spinon and a holon recombine into an electron 
at low energies 0, |2^|. Due to the spin-charge recom- 
bination, the pseudogap metallic state is described by 
electron-like quasiparticles at low energy. Since the bind- 
ing between the spinon and the holon is weak, the large 
pseudogap near the anti-nodal points (n, 0) and (0, n) 
is not affected. So the Fermi surface of the recombincd 
electrons cannot form a large closed loop. A simple the- 
oretical calculation^ suggests that the Fermi surface of 
the recombined electrons forms four small arcs near the 
nodal points (±tt/2, ±tt/2). Thus the SU(2) slave bo- 
son theory^ contains two key features: the pseudogap 




FIG. 6: Doping dependence of the truly measured su- 
perconducting transition temperature T c (full squares) and 
that calculated by T c = /3i>A7n(0) (open squares) with = 
0.7445-ft' 3 mo/s/ Jm. The solid line represents the empirical 
relation T c /T™ ax = 1 - 82.6(p - 0.16) 2 with T™ ax = 38K. 



due to spin singlet pairing and the Fermi arcs due to the 
spin-charge recombination[28j . And the superconductiv- 
ity arises from the coherent motion of the quasiparticles 
on the arcs, thus one expects that T c is proportional to 
the gap on the Fermi arc: k B T c « ^VAfok arc , instead 
of the pseudogap A p near the anti-nodal points. Mean- 
while, since the spin pairing is responsible for both the 
pseudogap A p near the anti-nodal points and the nodal 
gap slope VAi it is reasonable to see the proportionality 
between va and T*(cx A p ) or A q w 0A6k B T*. These are 
exactly what we found in the experiment. 



CONCLUDING REMARKS 

In summary, the Volovik's relation of the d-wave 
pairing symmetry has been well demonstrated by low 
temperature specific heat in wide doping regime in 
Lai- x Sr x CuO^. Based on this analysis the nodal gap 
slope va is derived and is found to follow the same dop- 
ing dependence of the pseudogap A p . This strongly indi- 
cates the close relationship between the pseudogap and 
the superconductivity. Meanwhile it is found that the su- 
perconducting transition temperature T c is determined 
by «A7n(0) instead of va- This discovery may suggest 
the importance of Fermi arcs near the nodal region and 
the superconductivity is induced by the formation of a 
new gap on these arcs. Both observations are consistent 
with the SU(2) slave boson theory based on the general 
RVB picture. 
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